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Abstract 

Analysis of data from the EGRET 7-ray detector on the Compton Gamma Ray 
Observatory indicated an anomaly in the form of an excess diffuse galactic flux 
at GeV energies over that which was theoretically predicted. Various explanations 
for this anomaly have been put forth, including the invocation of supersymmetric 
dark matter annihilation. We reexamine these explanations here, including a new 
discussion of the possible systematic errors in the sensitivity determination of the 
EGRET detector. We conclude that the most likely explanation of the EGRET 
"GeV anomaly" was an error in the estimation of the EGRET sensitivity at energies 
above ~1 GeV. We give reasons why such a situation could have occured. We find 
evidence from our new all-sky analysis which is inconsistent with the assumption 
that the anomaly can be a signal of supersymmetric dark matter annihilation. We 
also reconfirm the original results of the EGRET team on the extragalactic 7-ray 
background spectrum. There are important implications of our analysis for the 
upcoming Gamma Ray Large Area Telescope (GLAST) mission. 
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1 Introduction 



The EGRET 7-ray detector, a spark chamber telescope flown aboard the 
Compton Gamma Ray Observatory satellite, provided detailed measurements 
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of astrophysical 7-ray sources and galactic and extragalactic diffuse 7-ray 
fluxes. The reported galactic diffuse flux was measured and mapped over the 
whole sky. 

Theoretical studies of the physics and astrophysics of galactic 7-ray production 
provided predictions of the expected fluxes and energy spectra [1]. In the 
energy range above ~1 GeV , the fluxes reported by the EGRET team were up 
to ~60% higher than the theoretical predictions [2]. This apparent discrepancy 
between the reported fluxes and the theoretical calculations will be referred 
to here as the "GeV anomaly" . 

Unresolved galactic point sources can be ruled out as an explanation for the 
GeV anomaly for two reasons: (f) they would be strongly concentrated in the 
galactic plane and the GE anomaly, as we will show, is seen isotropically over 
the whole sky, and (2) the largest class of galactic point sources are pulsars 
and such sources are concentrated in the inner galaxy and they make up less 
than 15% of the total galactic flux [3]. 

There have been three approaches for accounting for the GeV anomaly as a 
diffuse phenomenon, viz. either (1) invoking a cosmic-ray electron source spec- 
trum proportional to E~ 2 with a resulting significant Compton 7-ray com- 
ponent above 1 GeV increasing the flux [4], (2) making modifications to the 
assumptions regarding both the primary cosmic-ray nucleon and electron spec- 
tra in numerical models in order to push up the total theoretical 7-ray fluxes 
[5], or (3) postulating new physics, namely supersymmetric dark matter an- 
nihilation, to account for the anomaly [6]. 

The first approach postulates that the cosmic-ray electron spectrum observed 
at Earth is much steeper than the average spectrum in the Galaxy, but that 
this ad hoc situation can result as an effect of the distribution of the supernova 
remnants (SNR) which produce the electrons and electron propagation effects. 
A prediction of this model is a center-anticenter asymmetry in the anomaly 
owing to the strong galactocentric radial distribution of SNR. A reduced bump 
at high galactic latitudes would also be expected, owing to the steeper local 
cosmic-ray electron spectrum. 

Approaches (2) and (3) have significant implications. Increasing the galactic 
7-ray production rate at GeV energies in the model of Ref. [5] produced a re- 
duction in the implied extragalactic diffuse flux in the followup calculation [7], 
with a marked dip at energies near ~ 1 GeV. This is opposed to the original 
determination of the extragalactic background spectrum by the EGRET team 
[8]. Postulating that the GeV anomaly is caused by a 7-ray component from 
supersymmetric dark matter annihilation has, of course, much greater implica- 
tions. It is therefore important to reexamine the issue of the GeV anomaly by 
paying careful attention to the collateral implications of deviations from the 
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canonical predictions [1] and to examine the more prosaic possibility that the 
assumed high-energy sensitivity of the EGRET detector may not have been 
correct. In this paper, we will make the case that a problem in the analysis 
of the EGRET sensitivity calibration is the most likely explanation for the 
"GeV anomaly" and that the dark matter hypothesis can be ruled out by an 
examination of the GeV spectrum over the whole sky. 



2 Production of Diffuse Galactic 7-rays 

Galactic 7-rays are produced by interactions of relativistic electrons and pro- 
tons with interstellar gas and photons. The physical processes involved are 
electron bremsstrahlung, neutral pion production, and Compton interactions 
of cosmic ray electrons [9] , [10] , [1] . Of these processes, the decay of neutral 
pions produced by cosmic rays interacting with interstellar gas is expected to 
dominate at energies above 0.1 GeV [1],[2]. The vast majority of these 7-rays 
are produced by cosmic rays with energies below ~ 20 GeV [11]. In this energy 
range the cosmic ray spectrum, particularly in the local galactic neighborhood 
that accounts for the 7-ray production at high galactic latitudes, is well mea- 
sured [12]. The pion production cross section at these energies is also very well 
known [11], [13]. At 7-ray energies up to ~1 GeV, the pion production process 
is well described by the "isobar-plus-fireball" model [9], [14], [15]. At energies 
above ~1 GeV, it is generally assumed that scaling holds and that the 7-ray 
spectrum will have the same spectral index as the primary proton spectrum 
[15], [16]. More recent calculations [18] are in good agreement with the results 
presented in Refs. [15] and [16], confirming that the physics of pion production 
at GeV energies is well determined. A modification of this physics involving 
scaling violation has been suggested [17], however the resulting effect is too 
small to explain the GeV anomaly by scaling violation alone. 

It has been suggested that the GeV anomaly can be explained by postulating 
that the average primary proton spectrum in the galaxy is harder than that 
observed locally. With scaling assumed, the required proton spectral index 
is T = 2.45 [18], [19]. Even with some scaling violation [17], one requires a 
spectral index T ~ 2.5. The problem with this assumption is that the local 
proton spectrum has a significantly steeper measured index T = 2.76 ± 0.03 
[12]. At galactic latitudes above ~ 20° from the plane, the diffuse 7-rays 
are produced by these locally measured protons (and a smaller number of 
cosmic-ray a particles which have a similar spectral index [12]). Even should 
the proton spectrum in the inner galaxy be harder than the local spectrum 
at high galactic latitudes, it is the steeper local spectrum which produces the 
high-latitude pion-decay 7-rays. Therefore, the predicted spectrum should be 
the canonical one. In other words, there should be no GeV anomaly at high 
galactic latitudes. This conclusion would also result if the anomaly is strictly 
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produced by a harder cosmic-ray electron spectrum averaged over the whole 
galactic disk [4]. 

On the contrary, it is generally accepted that the GeV anomaly is uniform 
over the entire sky, a result which was obtained from the EGRET data and 
which was implicit in Figs. 3 and 5 of Ref. [8] by the EGRET group. We will 
show here, in a very quantitative way, that the anomaly is seen at all galactic 
latitudes independent of the line-of-sight density of galactic hydrogen gas. We 
will further argue that such complete isotropy is what would be expected if 
this anomaly is most likely traceable to the detector itself. 



3 The Celestial Distribution of the GeV Anomaly 

Our detailed analysis of the galactic 7-ray flux has now been expanded to 
cover the entire sky [22]. This was accomplished by extending the latitude 
range of the galactic plane analysis to ± 25°, chosen to match the latitude 
range which includes all significant emission by molecular hydrogen clouds [21], 
and further extending the model to the galactic poles [22] by using the all- 
sky Leiden-Dwingaloo [23] and Instituto Argentino de Radioastronoma (IAR) 
Southern Hemisphere HI survey [24]. The three dimensional cosmic-ray den- 
sity for the latitude range \b\ < 25° was derived from the Galactic plane matter 
distribution on the assumption of dynamic balance [25]. The 7-ray production 
function per H atom in the 22° < \b\ < 25° latitude range was extended to the 
poles on the assumption that 7-ray production in this range is dominated by 
interactions with cosmic rays having the locally measured spectrum. The anal- 
ysis approach used in Ref. [8] to determine the extragalactic diffuse emission 
was then repeated using the new all-sky model for each of the ten standard 
energy ranges as well as the four broad energy ranges discussed in Ref. [2]. 
This analysis, using the EGRET phase 1 and phase 2 source-subtracted data, 
yields an extragalactic 7-ray spectrum which is consistent with the one derived 
in Ref. [8]. This consistency is expected because the extragalactic flux is the 
residual observed flux as the galactic model flux is extrapolated to zero. The 
slope of the extrapolation does not affect the intercept, as shown in Figure 1. 
Furthermore, this all-sky analysis confirms that the GeV anomaly is uniform 
over the entire sky, a result which was implicit in Figs. 3 and 5 Ref. [8], but 
which is now more quantitatively shown in Fig. 1. The lack of any structure in 
the anomaly related to the Galactic plane, galactic center, anti-center or halo, 
strongly indicates that the GeV anomaly is due to a systematic error in the 
EGRET calibration (see below) rather than being a real astronomical effect. 
The excellent 1:1 correlation of the EGRET data with the expanded all-sky 
model over the entire sky and the entire range of emission, after multiplica- 
tion by a single renormalization factor, further corroborates this conclusion 
as also shown in Fig. I. Table 1 and Fig. 2 breaks down this effect as a func- 
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tion of energy which indicates a miscalibration which gets worse at the higher 
energies. 

The overall effect of applying this renormalization correction for point sources 
is slight over a large energy range on a log-log plot and does not conflict with 
any theoretical models within the large observational error bars of the data 
above 1 GeV. For 7-ray blazars observed by EGRET, the fit to power-law 
spectra would result in slightly steeper spectra, with the error bars on specific 
points generally large in the energy range above 1 GeV. 

To specifically illustrate the effect of our renormalization factors on the spectra 
of point sources, we have applied our renormalization to the EGRET data on 
both the pulsed and unpulsed nebula spectrum of the Crab Nebula as shown 
in Figure 3 based on the compilation in Ref . [26] . The renormalized spectrum 
is consistent with previous theoretical models of 7-ray emission from the Crab 
Nebula [26], [27]. 



4 The Dark Matter Hypothesis 

It has been suggested that the GeV anomaly is the result of a new component 
to the diffuse 7-ray spectrum, one that is in addition to the pion-decay compo- 
nent, namely a component from supersymmetric dark matter annihilation [6]. 
However, there are serious problems with this hypothesis: (1) The same process 
of dark matter annihilation would produce a flux of cosmic-ray antiprotons 
which is incompatible with the measured value [28]. (2) The celestial 7-ray 
distribution from dark matter annihilation would be both highly asymmetric 
and clumped [29]. This is again in direct contradiction with isotropy of the 
anomaly. 



5 The Energy Calibration of the EGRET Detector 

It was pointed out in the EGRET calibration paper [30] that at high energies 
the EGRET sensitivity is poorly known, partly due to the uncertain effect 
of self-vetoing caused by charged-particle backsplash impinging onto the an- 
ticoincidence scintillator following shower generation in the detector [8]. This 
effect eliminates good 7-ray events. Any uncertainty in the self- veto correction 
would alter the sensitivity of the detector at energies above ~ 1 GeV. 

An attempt to quantify the amount of backsplash was performed experimen- 
tally up to 10 GeV during the calibration of EGRET at the Stanford Linear 
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Accelerator Center (SLAC) with over a thousand hours of 7-ray beam ex- 
posure. The beam was scanned across the face of the tracking detector at a 
series of off-axis angles from 0° to 40° with three different roll-axis angles and 
10 different energies. An EGS (Electron Gamma Shower) Monte Carlo (MC) 
code was used to study the effect of the backsplash under differing conditions, 
but the fidelity of the MC model was low. It did indicate that the backsplash 
was strongly dependent on the charged-particle threshold, which was known to 
be strongly variable across the anti-coincidence detector. The photomultiplier 
gains were reduced after the SLAC calibration to mitigate the higher than 
anticipated backsplash effect. The effect of lowering the gains is substantial, 
and is not reflected in the sensitivity calibration results as given in Ref. [30]. 

The all-sky pervasiveness of the apparent excess flux above 1 GeV energy 
(see above) strongly indicates that this anomaly is intrinsic to the detector. A 
systematic error in the understanding of the detector response would be one 
likely source. In an effort to investigate this possibility, we have reexamined 
the EGRET archival data carefully, selecting only the highest quality EGRET 
data. The selection was based on: (1) choosing 7-ray events with the most reli- 
able energy determination, (2) selecting 7-ray events within 30° of the detector 
axis for which the instrument response is most reliably determined, and (3) ex- 
cluding 7-ray events within a 4cr angular uncertainty from the Earth's horizon. 
We have examined data on the pulsed and unpulsed emission from the Crab, 
Vela, and Geminga and the diffuse flux from the Lockman Hole for phases 
1 - 4 of the EGRET data for energies > lGeV, We find that for different 
observations of these sources, there were variations in their measured fluxes 
in excess of 40%. Since the flux from these sources is expected to be constant, 
this indicates systematic errors caused apparent time variations over different 
observing periods. This is not surprising, since the effect of degraded perfor- 
mance with the aging of the spark chamber gas was calibrated in flight by 
fixing the > 100 MeV flux and not taking account of the energy dependence 
of this effect [31], thus producing a false variability at GeV energies. 

While we find anomalous behavior pointing to an imprecise knowledge of the 
detector response, we are unable to explain the GeV anomaly as caused by any 
single systematic instrument effect that we have studied. Other contributors 
besides an imprecise knowledge of the detector response that could explain 
an apparent trans-GeV anomaly include (1) an error in the analysis which 
might cause non- 7-ray events to be identified as good high-energy events, (2) 
an incorrect assignment of energies and (3) an induced detector background 
in the harsh environment of a low-Earth orbit. It is difficult to eliminate any 
of these possibilities, but no specific evidence has been found that they exist. 

Given the systematic energy dependent uncertainties in the sensitivity cal- 
ibration for the EGRET instrument discussed in this paper, the published 
EGRET results above 1 GeV cannot be reliably depended upon. We suggest 
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Table 1 

INVERSE NORMALIZATION FACTOR VS. ENERGY 



Energy Bin 


Inverse Normalization Factor 


30-50 MeV 


0.815 ± 0.009 


50-70 MeV 


0.992 ± 0.018 


70-100 MeV 


0.947 ± 0.019 


100-150 MeV 


0.959 ± 0.019 


150-300 MeV 


0.925 ± 0.013 


300-500 MeV 


0.887 ± 0.022 


500-1000 MeV 


0.965 ± 0.026 


1-2 GeV 


1.26 ± 0.06 


2-4 GeV 


1.57 ± 0.17 


4-10 GeV 


1.85 ± 0.46 



that one should incorporate the energy dependent correction factor shown in 
Figure 2. This suggestion applies to possible use by the GLAST LAT (Large 
Area Telescope) collaboration for comparison as a sensitivity calibration check 
on their extensive Monte Carlo simulations and beam tests. We must await 
results from GLAST in order to pursue this question further. 



6 Conclusions 



Our new all-sky analysis of the EGRET data confirms a systematic anomaly 
in the form of an apparent excess galactic emission at energies above 1 GeV 
over what is expected from galactic cosmic ray interactions with gas nuclei 
in atomic and molecular clouds and interstellar photons, as previously re- 
ported [2]. We find further quantitative support that this anomaly is constant 
over the whole sky, not being correlated with any astronomical or galactic 
features. It is therefore most likely caused by a systematic error in the cali- 
bration of the effective detector sensitivity. Although a detailed reanalysis of 
the calibration data is impossible at this time, we have shown that plausible 
systematic uncertainties in the calibration of the EGRET sensitivity for 7-ray 
energies above 1 GeV can readily account for the universal anomalous excess 
flux. We conclude that neither making modifications to the observed primary 
cosmic-ray nucleon and electron spectra in order to push up the predicted 
theoretical 7-ray flux [5] nor postulating a component from supersymmetric 
dark matter annihilation [6] are required to account for the "GeV anomaly" . 
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Our new analysis is consistent with the extragalactic flux derived in Ref. [8]. 
This reevaluation of the EGRET sensitivity also has implications for calibra- 
tion checks of the GLAST 7-ray detector to be launched in the near future. 
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Fig. 1. (a) Plot of integral (E > 1 GeV), all-sky diffuse model flux vs. EGRET 
observed flux for 335° < I < 45° , |6| < 90°. (b) A similar plot with a renormal- 
ization factor of (1.6) -1 applied to the observed flux. In both plots, the integral 
extragalactic diffuse flux of 1.5 xlO -6 cm -2 s _1 sr _1 has been added to the diffuse 
model. The dotted line indicates the expected 1:1 relationship between the model 
and observed fluxes. Contours show the number of 0.5° x 0.5° pixels containing 
the flux indicated within a bin of width 1.6 x 10~ 6 cm -2 s -1 sr _1 . The contour val- 
ues are 10 4 , 10 3 , 10 2 , 10 1 , ... These plots clearly demonstrate that the GeV anomaly 
exists uniformly over the whole sky and extends from high to low intensity galactic 
flux emission. 
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Fig. 2. Required inverse renormalization factor for different energy bins, given as 
the ratio of observed-to-predicted flux vs. energy. 
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Fig. 3. Pulsed and Nebula spectrum of the Crab Nebula with and without the cor- 
rection factors from Figure 2 applied to the EGRET data, based on the compilation 
in Ref. [26]. The crosses in the energy range from 50 MeV to 10 GeV show the 
original EGRET data points and the diamonds show the renormalized values. 
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